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ABSTRACT 

We show in this paper why molecular millimeter absorption line searches in DLAs have been unsuccessful. 
We use CO emission line maps of local galaxies to derive the H2 column density distribution function f(Nn 2 ) 
at z = 0. We show that it forms a natural extension to f(Nm): the H2 distribution exceeds /(iVm) at 
Njj ~ 10 22 cm~ 2 and exhibits a power law drop-off with slope ~ —2.5. Approximately 97% of the H2 mass 
density pn 2 lS m systems above 7Vh 2 = 10 21 cm" 2 . We derive a value pn 2 = 1-1 x 10 7 /i7o M Q Mpc~ 3 , which 
is « 25% the mass density of atomic hydrogen. Yet, the redshift number density of H2 above this Nn 2 limit 
is only ss 3 x 10~ 4 , a factor 150 lower than that for H I in DLAs at z = 0. Furthermore, we show that the 
median impact parameter between a 7Vh 2 > 10 21 cm" 2 absorber and the centre of the galaxy hosting the H2 
gas is only 2.5 kpc. Based on arguments related to the Schmidt law, we argue that H2 gas above this column 
density limit is associated with a large fraction of the integral star formation rate density. Even allowing for an 
increased molecular mass density at higher redshifts, the derived cross-sections indicate that it is very unlikely 
to identify the bulk of the molecular gas in present quasar absorption lines samples. We discuss the prospects 
for identifying this molecular mass in future surveys. 

Subject headings: galaxies: ISM — ISM: molecules — ISM: atoms — quasars: absorption lines 



1. INTRODUCTION 

At high redshift, the atomic hydrogen is well surveyed 
through observations of damped Lya systems (DLAs), quasar 
absorption line systems characterized by -/Vhi > 2 x 
10 20 cm -2 (e.g.. lProchaska et al.l l2005). These observations 
yield the frequency distribution of H I surface density /(A/hi) 
and the first moment gives the cosmological mass density 
of predominantly neutral, atomic gas. Several small surveys 
have been performed to search for molecular gas associated 
with DLAs. Yet, there have been no positive detections of 
CO and other m olecules in m illimeter wavelength absorp- 
tion line searches (Curran et al. 2004, and references therein). 
Searches for H2 absorption in DLAs (via the Lyman and 
Werner bands) show a success rate of only w 20% and even 
these sightlines have low molecular fractions (Ledo uxet al.l 
2003). 

These low detection rates have drawn into question the 
relationship between DLAs an d star formation at high red- 
shift (e.g. Honkins et al. 2005). Other galaxy surveys have 
shown that star for mation is very active at these epochs (e.g., 
Steide fet alJl99 9). Presumably, since stars form in molecular 
clouds, these star forming galaxies have significant molecular 
gas mass that has not been discovered through quasar absorp- 
tion line studies. The obvious conclusion, then, is that these 
molecular regions have very low covering fraction on the sky 
and/or contain sufficient columns of dust to obscure any back- 
ground quasar. 

In this paper we address these issues by makin g use of the 
studies of the CO distribution in nearby galaxies. Zwa an et al.l 
( 2005b) recently showed that most DLA properties (luminosi- 
ties, impact parameters between quasars and DLAs, and metal 
abundances) are consistent with them arising in galaxies like 
those in the local universe. Building on this result, we use 
here information on the molecular content of nearby galax- 
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ies, to make predictions of the detectability of H2 in DLAs. 
The H2 distribution in nearby galaxies is generally derived 
through mapping of the CO(1-0) line at a frequency of 115 
GHz. The H2 column densities derived from these observa- 
tions are typically higher than 5 x 10 20 cm -2 , which is very 
close to those reached by millimeter CO absorption line ob- 
servations in redshifted DLAs against radio-loud background 
sources. 

2. THE H 2 CROSS-SECTION 

The largest sample of high spatial resolution CO(1-0) maps 
of nearby galaxies is the BIMA SONG sample presented by 
Heifer et al. (2003). The sample consists of 44 optically se- 
lected galaxies with Hubble types Sab to Sd at a median dis- 
tance of 11.9 Mpc, and is selected without reference to CO 
or infrared brightness. The fact that this sample does not in- 
clude early types and irregulars is not expected to introduce 
an important bias in our calculations. CO has been detected 
in a number of elliptical galaxies ( Young|2002|), but their con- 
tribution to the total H2 mass density is expected to be mini- 
mal. Irregular s contr ibute approximately 20% to the total H I 
mass density (|Zwaan et alj l2003l) . but their ratio of H2 mass 
to H I mass is down by almo st a factor ten compare d to the 
average over all galaxy types (lYoung & Knezekfl9 89). which 
implies that their contribution to t he H2 mass density is also 
marginal. We refer to Zwaan et al. ( 2005b) for a discussion on 
why optically selected galaxies are a fair sample to compare 
with absorption line statistics. 

We used the combined interferometric and single dish maps 
to calculate H2 column density maps, using a CO/H2 con- 
versio n factor of 2.8 x 10 20 cm~ 2 (Kkms" 1 ) 1 l lKennicuttl 
1998). iHelfer et alJ J2003I) report that the lowest signifi- 
cant isolated H2 densities at 3<r in these maps are approxi- 
mately 13.7 Mq pc -2 , which corresponds to an 7Yh 2 limit of 
8.5 x 10 20 cm~ 2 . The typical spatial resolution of the 6" beam 
is 350 pc at the median distance of the sample. 

We calculate the H2 column density distribution f(Nn 2 ) 
using a method analogous to that for Hi set out by 
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IZwaanetalJll2005bl) : 

f(N x_ c Y.iHMB l )w(M Bl )A. l {\ogN Yl2 ) 
J{ Ha) H N H2 lnlOAlogiVH 2 ' ^ ' 

Here, $(Mbj) is the space density of galaxy i measured 
through the optical luminosity function as measured by 
Norb erg et alJ (120021) . with Schechter parameters: M* B — 
51og/i 70 = -20.43, a = -1.21, and $* = 5.5 x 
10 _3 /i 3 Mpc~ 3 . The function w(MBi) is a weighting func- 
tion that takes into account the varying number of galaxies 
across the full stretch of Mb, and is calculated by taking the 
reciprocal of the number of galaxies in the range Ms, — A/2 
to M Bi + A/2, where A is taken to be 0.25. Aj(logiVH 2 ) is 
the area function that describes for galaxy i the area in Mpc 2 
corresponding to a column density in the range logAn 2 to 
log N-r 2 + A log N-r 2 . In practice, this is simply calculated 
by summing for each galaxy the number of pixels in a certain 
log N-r 2 range multiplied by the physical area of a pixel. Fi- 
nally, c/Hq converts the number of systems per Mpc to that 
per unit redshift. 

The BIMA SONG galaxies are selected to be less inclined 
than 70°. In order to achieve a f(N) measurement for ran- 
domly oriented galaxies, we de-projected all galaxies to face- 
on assuming that the H2 gas is optically thin, and subse- 
quently recalculated the column density distribution function 
for ten inclinations i evenly spaced in cos(i) between and 1. 
The final f(Nn) was taken to be the average of these ten dis- 
tribution functions. This procedure only makes a small mod- 
ification to the f(Nji) calculated from the H2 distributions 
uncorrected for inclination effects. 

Figure [2 shows the resulting column density distribution 
function f(Nn 2 ), together with f(Nm) from iZwaan et all 
(2005b). We note that the horizontal axis represents the atom 
surface density, which in the case of H2 is equal to 2Ah 2 . The 
errorbars are la uncertainties and include counting statistics 
and the uncertainty in the optical luminosity function. The un- 
certainty in the CO/H2 conversion factor could introduce the 
largest error in our /(Ah 2 ). The horizontal errorbar indicates 
the uncertainty in f(Nn 2 ) if this conversion factor is uncer- 
tain by 50%. The f(Nn 2 ) can be fitted very well with a log- 
normal distribution: f(N) — f* exp -[(log N — a)//j] 2 /2 
, where /1 = 20.6, a — 0.65 and the normalization /* is 
1.1 x 10 _25 cm 2 . The distribution function of H2 column den- 
sities seems to follow a natural extension of the H I distribu- 
tion function, in such a way that the summed f(Nn) roughly 
follows a power-law distribution N^ 2 ' 5 between log Ah = 
21 and log Ah = 24. The two distribution functions cross at 
log Ah ~ 22, which is the approximate column densit y asso - 
ciated with the conversion from H I to H2 (e.g. Schave 2001 ). 

By taking the integral over /(Ah 2 ) multiplied by Ah 2 , 
we find the total H 2 mass density at z = to be pn 2 
1.1/ 10 7 M (7i Mpc~ 3 , which is approximately one qu arter of 
p m (z = 0) llZwaan etai1l2005al) . iKeres et all (120031) found 
pn 2 = (2.0 ± 0.7) x 10 7 M Q Mpc -3 , where we converted 
their value to Ho = 70 km s -1 Mpc -1 and used the same 
CO/H2 conversion factor as we used. This difference might 
be partly due to the fact that the BIMA SONG sample does 
not include many d warf galaxies, although judging from the 
IKeres etall ll2003l) CO mass function, these galaxies con- 
tribute only very little to pn 2 • Another reason could be that 
the BIMA SO NG sample is optica lly selected, whereas Keres' 
sample (Young & Knezekl ll989l) predominantly consists of 
FIR-selected galaxies, which can cause a bia s toward CO- 
rich galaxies (see e.g.. lSolomon & Sagelll988l) . ICasoli et all 
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FIG . 1 . — The column density distribution function of H I and H2 at z = 0. 
The Hi curve is adopted from Zwaan et al. 1 2005b), the H2 curve is mea- 
sured from the CO emission line maps from the BIMA SONG sample. Col- 
umn densities are expressed in atoms per cm , also for H2 . The solid line 
is the summed /(-/Vjj)- The horizontal errorbar indicates the uncertainty in 
f(Nft 2 ) if the CO/Nft., conversion factor is uncertain by 50%. The /(iVu, ) 
can be fitted very well with a log-normal distribution, where /i = 20.6, 
a = 0.65 and the normalization is 1.1 X 10 — 25 cm 2 , as indicated by the 
dashed line. 



(1998) used a larger sample than Youn g & Knezekl (Q989), 
and took into account CO non-detections, to find much lower 
values of Mr 2 /Mhi, which are fully consistent with our de- 
rived value of ph 2 / Phi- 

What fraction of the cosmic H2 mass density do we miss be- 
low the A H2 detection limit of 8.5 x 10 20 c m - 2 ? The /(A H J 
appears to flatten off at the lowest column densities, which 
implies that the contribution of low Ah 2 is low, unless the 
f(Nn 2 ) rises steeply below our detection limit. To test this, 
we make use of the H2 absorption line survey in DLAs by 
Ledoux et al. (2003). These authors report a detection rate of 
13 to 20%, with H2 column densities typically in the range 
log Ah 2 = 17 to 18.5. Based on their statist ics and the mea - 
sured /(A H i) of DLAs at z = from lZwaan et alJ J2005b). 
we estimate that log /(Ah 2 = 10 18 cm~ 2 ) = -23.5 at z = 0, 
conservatively assuming that the detection statistics of H2 ab- 
sorption has not evolved since z « 2. Extrapolating our mea- 
sured /(AW) to this value, we find that /(A H2 ) oc A H2 ~ 0,5 
(see Figure|2). Now, by integrating f Ah 2 /(Ah 2 )c?A 7 h 2 we 
find the total H2 mass as a function of Ah 2 . From this we 
derive that the mass contained in systems logAn 2 < 21 is 
only 3% of the total H2 mass density. This implies that our 
results presented in this paper apply to roughly 97% of the 
total number of H2 molecules in the universe. 

Beam smearing might lead to an overestimation of the 
cross-sections at low Ajj 2 and an underestimte of f(Nn 2 ) at 
large Ah 2 . Because a main result of the next section is that 
the H2 cross-sections are small, we have ignored these effects 
for now. 

3. WHERE TO FIND THE MOLECULES 

The redshift number density dN/dz of H2 above the Ah 2 
limit of 10 21 cm~ 2 can be calculated from f(Nn 2 ) by sum- 
ming over all column densities larger than this limit. We find 
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FIG. 2. — Top: The column density distribution function of H2, where the 
circles show f(Nn 2 ) at 2 = derived from the BIMA SONG CO maps, and 
the grey square indicates the value at much lower column densities, estimated 
from the H2 absorption data from Ledoux et al. 1 2003). These latter data are 
from systems at redshifts 2^2, which implies that the derived f(Nn 2 ) 
point is probably an upper limit. Bottom: The cumulative H2 mass distri- 
bution in column densities > -/Vjj 2 ■ Here, we have includes the -/Vh 2 — °' 5 
extension as shown in the top panel. Approximately 97% of the H2 molecules 
are in column densities Afjj 2 > 10 21 cm -2 . 



that dN/dz = 3 x 10~ 4 , which is approximately 150 times 
lower than the corresponding value for H I above the DLA 
limit fZwaan et al. 2005b). Taking this result at face value, 
this would imply that at z = only one in every ^150 DLAs 
is expected to show strong CO absorption lines. 

In Figure [3] we show how close to the centre of galaxies 
high column density H2 is to be found. This figure shows 
the normalized cumulative distribution function of impact pa- 
rameters between the centers of galaxies and the positions 
where the H2 absorption is measured. We constructed this di- 
agram from the BIMA SONG CO maps, and using a weight- 
ing scheme similar as discussed for the evaluation of f(Nn 2 ). 
Apparently, the median impact parameter at which an 7Yh 2 > 
10 21 cm~ 2 absorber is expected is only 2.5 kpc and 90% of 
the cross-section is at impact parameters smaller than 6.5 kpc. 
For higher H2 column densities, these impact parameters are 
even smaller. Applying these findings to the DLA popula- 
tion, millimeter CO observers need to select those DLAs that 
arise within 2.5 kpc of the centers of galaxies to have a 50% 
probability of identifying an absorption line corresponding to 
Nn 2 > 10 21 cm~ 2 . None of the identified DLA host galaxies 
have such small impact parameters, except perhaps two cases 
studied by Rao et al. ( 2003) for which only upper limits to the 
impact parameters could be given due to blending with the 
background quasar light. 

Our conclusions are qualitatively supported by observations 
of redshifted molecular lines. To date, only four reds hifted 
molec ular absorbers have been identified ( Wiklind & Combes 
1999, and references therein), two of which arise in gravita- 
tionally lensed intervening galaxies and two originate in the 
galaxy that also hosts the radio source against which the ab- 
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FIG. 3. — The normalized 2 = cumulative distribution function of impact 
parameters of H2 cross-section selected systems above different H2 column 
density limits. The solid is the distribution for impact parameters < b, the 
dotted line for > b. 



sorption is seen. In all cases the sightline through the high H2 
column passes through the galaxies at small impact parame- 
ters. 

These cross-section arguments already explain the very low 
detection rate of molecules in DLAs, but there are additional 
effects that hinder the detection of molecules. Molecular hy- 
drogen forms on the surface of dust grains. The regions in 
galaxies containing most of the universe's H2 molecules are 
therefore likely to be dusty, causing a higher optical extinc- 
tion of the background sources, which, in turn, might lead to 
them dropping out of magnitude-limited surveys. Obviously, 
radio-selected quasars do not suffer from this bias and are po- 
tentially good candidates against which to find high iVH 2 ab- 
sorbers. However, dN/dz(Nn 2 ) is so low that identifying 
such an absorber is highly unlikely in current radio-selected 
quasars samples. For exa mple, the reds hift interval covered 
by the CORALS survey JEllison etalJl2001l) of radio-loud 
quasars is only Az = 55 for Lya, corresponding to Az w 100 
for H2. The redshift interval required for identifying a high 
N-r 2 absorber should be on the order of (dN/dz) -1 — 3300 
at z = 0, or approximately 600 at z — 2 — 3, taking into 
account the cosmological variation of the absorption distance 
interval dX. Using this method to obtain sufficient statistics 
to measure the molecular mass density at intermediate and 
high redshifts, would require radio source samples larger than 
currently available. 

Finally, we address the question whether there might ex- 
ist a significant amount of H2 gas not associated with DLAs. 
CO mapping of nearby galaxies has shown that many galaxies 
show a depression in the H I distributi on where the H 2 col- 
umn densities are highest (see e.g. Wong & Blitz 2002). Pre- 
sumably, in regions where the molecular densities are highest, 
most of the H I has been converted to H2. We use the H I and 
CO m aps of seven nearby galaxies studied by Wons & Blitzl 
d2002l) to test how much H2 exists below the H I DLA limit. 
We compare the H I and CO maps smoothed to the same spa- 
tial resolution of 13" to 23", and find that the mass fraction 
of H2 in regions where iVni < 2 x 10 20 cm~ 2 ranges from 
to 40%. Averaged over all pixels from the total sample, this 
fraction is 6%. In many cases the fraction is underestimated 
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because the highest Ah 2 is smoothed over larger regions. Of 
course, this small sample may not be representative of the to- 
tal galaxy population at z = 0, but at least this exercise shows 
that a small fraction of the cosmic H2 density must be found 
in sub-DLAs. One might identify such absorbers by searching 
for a sub-DLAs with abnormally high metallicities, i.e. this 
could indicate absorbers where the hydrogen is not predomi- 
nantly atomic. An example of a high metall icity sub-DLA has 
recently been found by Perou xet alJ (12006). 

4. IMPLICATIONS FOR THE STAR FORMATION RATE DENSITY 



lLanzetta etaO (120021) and iHopkins etail (120051) recently 
estimated the star formation rate density (SFRD) in DLA 
systems by applying the 'Schmidt law' of star formation to 
the Hi column density distribution function /(Ahi). The 
Schmidt law is defined in local galaxies and states that the star 
formation rate correlates very well with total neutral gas sur- 
fac e density £m + ^ h, to the power 1 .4, as was demonstrated 
bv iKennicuttl ( 119981 . Here, we wish to investigate whether 
the SFRD measurements based on H I alone give meaning- 
ful results, or whether H2 should be taken into account for a 
reliable SFRD measurement. Since presently we only know 
f(N-Rn) at z = 0, w e cannot improve the measurements of 



Hopkins et al. (2005) at high redshift by including H2. Our 
aim is simply to test the validity of the method and discuss 
the implications of including molecules. 

We start with estimating what fraction of the SFRD is ac- 
tually contributed by those regions where the H2 column is 
higher than the Hi column. To this end, we make use of 
Equation 4 from Hopkins et al. (2005) that relates the SFRD 
to the f(Nu), an d apply it to our mea surement of /(Ah 2 ), 
and /(Ahi) from Zwaan et al. (2005b). We make one im- 
portant modification in that we convert our z — gas den- 
sities to those that would be observed if the gas disks were 
observed 'face-on', assuming that the Hi and H2 layers are 
optically thin. At z = the /(Ah) is a result of a randomly 
oriented population of galaxies and it is easy to see that the 
highest column densities are mostly the result of highly in- 
clined disks. However, the Schmidt law is valid for face-on 
gas densities. De-projecting the disk implies that the sur- 
face area increases but the column densities decrease. The 
net result is that that the /(Ah) as it is normally defined will 
overestimate the SFRD. At z = 0, we find that the SFRD 
is overestimated by ~ 40% if the regular /(-/Vh) is used, in- 
stead of the de-projected /(Ah) for face-on disks. The de- 
projected /(Ah) follows a nearly exponential behaviour be- 
tween log Ah = 20.5 and log Ah = 23.5, and can be fitted 
with the simple equation /(A H ) = 2 x 10 29 A^ 25 . 

Our results are presented in Fig. |4] which shows the im- 
plied SFRD as a function of H I and H2 face-on column 
density. We see that the H I and H2 column densities con- 
tribute approximately equally to the total SFRD. The other 
conclusion from this exercise is that the total derived SFRD 
at z = is much higher than that derived from Ha and 
rOul measure ments [p*(z = 0) rj 0.02M m yr" 1 Mpc~ 3 , 
Hopkins 2004]. This contradicts the findings of Hopkins et al. 
(2005) who conclude that at z = the SFRD p* derived 
from the Schmidt law and the measured /(Ahi) agrees well 
with the direct measurement s. The origin of this contra- 
diction is that Hopkins et al\ (|200 5|) use d the /(Ahi) mea- 
surements of iRTan^WebeTetalir iEoOS ). which were later 
corrected upwa rds with a factor 3. If we use the cor- 
rected values ( Rvan- Weber et al] 120051) or the measurement 
from IZwaan et alJ (l2005bl) . we find that p*{z = 0) = 
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FIG. 4. — The implied star formation rate de nsity as a functio n of face-on 
H I and H2 column density as derived from the Kennicutt 1 1998) star forma- 
tion law. Grey areas indicate approximate uncertainties. 



0.035M Q yr -1 Mpc~ 3 , a factor two higher than the me- 
dian of the directly measured values at z = 0. Us- 
ing the more realistic face-on /(Ah), and including both 
Hi and H2 in the analysis, we find that p*(z = 0) = 
0.044Mq yr -1 Mpc~ 3 , also higher than the nominal value 
at z = 0. Based on the face-on /(Ah 2 ) only, we find 
p*(z = 0) = 0.022Afoyr~ 1 Mp c~ 3 , in good agreeme nt with 
the direct measurements (see also Wong & Blitz 2002). 

Why is the value of p* at z = overestimated by a factor 
2 — 3 when derived using the /(Ah) and the Schmidt law? 
The answer lies in the fact that the Schmidt law is defined 
for the 'star-forming disks', and not for the region s of th e H I 
layer outside this area. It was shown by Kennicutt ( 1983) that 
star formation only occurs when the gas density exceeds the 
critical threshold density, which depends on the velocity dis- 
persion of the gas and the galaxy's rotation curve shape and 
amplitude. From the /(Ah) alone it is impossible to deter- 
mine what the critical threshold gas density is, because this 
density varies between galaxies and within galaxies. There- 
fore, applying the Schmidt law to all regions in the local uni- 
verse where the Ahi exceeds ~ 10 20 cm~ 2 will grossly over- 
estimate p*. [As an extreme example, consider NGC 2915, 
( Meur er ;t alJfl9 96) where the HI disk is many times larger 
than the optical disk.] 

Furthermore, within galaxies the areas with the highest 
SFRs are often mostly molecular, and in some cases Ahi 
declines in regions of high S FR (Martin & Kennicutt! 1200 It 
iRownd & Younj l Il999t |Wong & Blitzl 120021) . Inregions 
where log Ah 2 > 21 the critical density is typically exceeded. 
Consequently, for the purpose of estimating p*, the /(Ah 2 ) 
is probably a more reliable estimator, whereas /(Ahi) only 
gives an upper limit to that fraction of p* contributed by re- 
gions that are mostly atomic. 

In our analysis we treat the H I and H2 independently. Ide- 
ally, we would use H I and H2 measurements from the same 
large galaxy sample, but unfortunately such a sample is not 
available. Our analysis probably slightly wnt/erestimates the 
total SFRD: in regions where the Hi and H2 columns are 
equal, the SFR would be 2 14 /2 = 1.3 times higher if it were 
derived from the summed gas density instead of from the in- 
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dividual densities. In most regions, where the two column 
densities are different, the SFR is underestimated by a much 
smaller fraction if the H I and H2 are treated independently. 

In summary, we have shown here that /(iVni) cannot be 
used to derive a meaningful SFRD at z = and therefore 
probably also not at higher redshifts. Only in the unlikely situ- 
ation that DLA galaxies were truly molecule free, the f(Nm) 
method would give a meaningful upper limit to the SFRD con- 
tributed by DLAs. Interestingly, the sha pe of /(A/hi) has not 
evolved 3 between z = 4 and z = (cf. lZwaan et all 2005b: 
iProchaska et al.l2005l) . only the normalization (i.e., fim) has 
decreased wit h approximately a factor two. Based on the 
Honki nTet al.l d2005) parameterization, this implies that as 
measured through DLAs also has dropped only approximately 
a factor two over this redshift range. This is clearly at odds 
with the observations, unless one is willing to assume that the 
systems holding the bulk of the neutral gas account for only 
a small fraction of the cosmic SFRD. Furthermore, even at 
z = the SFRD is dominated by re gions where A^h 2 exceeds 
A/hi (see also [Wong & Blitz 2002). Given the observational 
fact that evolves much faster than SIhi, this implies that the 
cosmic H2 mass density was much higher in the past, or the 
laws of star formation were differ ent. 

In any case, the contention of Hop kins et al.l d2005l) that 
DLAs cannot make up for the observed SFRD at high red- 
shift should be reviewed. Clearly, to evaluate the SFRD con- 
tributed by DLAs one should take into account their molecu- 
lar content, which we have shown in the previous paragraph is 
extremely difficult to determine. In the absence of a measure- 
ment of f(Nn 2 ) at high er redshifts, more direct approaches 
such as the Cll* method (Wol fe et al.|2003l ) or models o f the 
star formation history (Dessauges-Zavadsk v et al.ll2003l) are 
probably more useful in constraining from DLAs. 

5. CONCLUSIONS 

We have used observations of CO in nearby galaxies to 
make predictions on the detection rate of H2 in DLAs. We de- 
rived a column density distribution function /(A/h 2 ) and find 
that it forms a natural extension of /(A?hi) at higher hydro- 
gen column densities. The inferred redshift number density 
dN/dz of A^h 2 > 10 21 cm" 2 is 3 x 10~ 4 , which is a factor 
150 lower than the corresponding number for DLAs. This im- 
plies that of the total number of DLAs currently known, only 
a handful are expected to show high column density molec- 
ular absorption lines. This absorption is expected to arise at 
impact parameters smaller than 6.5 kpc in 90% of the cases. 
Approximately 97% of the cosmic H2 mass density resides in 

3 Note however that Rao et al. 1 2005) find a much flatter f(Nm) distribu- 
tion at 0.1 < z < 1.7, based on Mgll-selected systems. 

4 http://www.alma.info/ 

5 http://www.skatelescope.org/ 



this dense, central gas, which also is associated with the bulk 
of the cosmic star formation rate density p*. The detectabil- 
ity of this molecular gas is further reduced by the facts that 
i) the high column density H2 is likely associated with high 
dust columns, which obscures background quasars and could 
take them out of magnitude-limited surveys, and ii) a small 
fraction (< 10%) of the H2 is associated with sub-DLA Lya 
absorption. 

We apply the Schmidt law for star formation to /( Nn 2 ) and 
/(A/hi) an d find that molecules should be taken into account 
to derive meaningful values for the cosmic star formation rate 
density p*. Even though the cross section of high column 
density H2 is low, this gas is associated with a large fraction of 
p». Since the molecular gas is difficult to detect in absorption, 
the method of using the Schmidt law and /(A/h) to derive 
remains of limited use. 

The cold atomic gas content of the universe at redshifts 
z ~ 1 — 5 has successfully been determined from blind sur- 
veys for DLA absorbers. What are the prospects for mea- 
suring the molecular mass density at intermediate and high 
redshifts using similar techniques? With present technology 
blind molecular absorption line surveys against radio-loud 
background sources seem difficult given the required redshift 
interval. The prospects for future instruments seem better. For 
example, with the Atacama Large Millimeter Array 4 (ALMA) 
the CO(3-2) line corresponding to log A/h 2 > 21 should be 
detectable at z = 2 to 3 in one minute against a 50 mJy 
background source. A survey of several days should be able 
to turn up a useful number of high A^h 2 absorbers. Such a 
survey would by simultaneously sensitive to HCO + (a good 
H2 tracer) at lower redshifts. The Square Kilometer Array 5 
(SKA) should be able to detect the CO(l-O) and HCO+(1-0) 
lines at redshifts z > 3.6 and > 2.6, respectively. The ex- 
pected noise levels are much lower than for ALMA, implying 
that this instrument will be ideal for high-z molecular absorp- 
tion line surveys, giving good statistics for measuring 17h 2 ■ 
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